Spin-orbit coupling and spectral function of interacting electrons in carbon nanotubes 
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The electronic spin-orbit coupling in carbon nanotubes is strongly enhanced by the curvature 
of the tube surface and has important effects on the single-particle spectrum. Here, we include 
the full spin-orbit interaction in the formulation of the effective low-energy theory for interacting 
electrons in metallic single-wall carbon nanotubes and study its consequences. The resulting theory 
is a four-channel Luttinger liquid, where spin and charge modes are mixed. We show that the 
analytic structure of the spectral function is strongly affected by this mixing, which can provide an 
experimental signature of the spin-orbit interaction. 

PACS numbers: 71.10.Pm, 72.15.Nj, 73.63.Fg, 72.25.Rb 



Spin-orbit interaction (SOI) effectsi are of great in- 
terest in the field of spintronics, and their detailed un- 
derstanding is both of fundamental and of technologi- 
cal interest, e.g., for the coherent manipulation of spin 
qubitsi^ In single-wall carbon nanotubes (SWNTs) the 
SOI arises predominantly from the interplay of atomic 
SO coupling and curvature-induced hybridization, and 
its effect on the electronic bandstructure has recently 
been clarifiedi^^— Contrary to other carbon-based ma- 
terials, as, e.g., flat graphene, where SOI is very weak 
(of the order of few /xeV), in SWNTs it can reach frac- 
tions of meV and has important consequences, previously 
overlooked. New experimentsi^ii^ on ultraclean SWNTs, 
made possible by advances in the fabrication technology, 
have indeed observed modifications of the electronic spec- 
trum due to SOI. These results confirm the theoretical 
expectations, and motivate renewed interest on SOI in 
nanotubes. So far on the theory side the main focus 
has been on nanotube quantum dots;^ where the SOI 
manifests itself in spectral features. Here we study long 
SWNTs, where long-ranged interactions can induce non- 
Fermi liquid electronic phases In particular, without 
SOI, Luttinger liquid (LL) behavior— has been predicted 
for metallic SWNTs Experimental evidence for this 
strongly correlated phase has been reported using quan- 
tum transporti^ and photoemission spectroscopy^^ 

The question therefore arises how the LL theory of 
SWNTsii is modified when the SOI is taken into ac- 
count, and what are its observable consequences. This 
question is addressed and answered in detail below. Our 
main results are as follows, (i) The low-energy theory of 
metallic SWNTs still describes a Luttinger liquid. How- 
ever, the decoupled plasmon modes do not correspond 
to spin and charge anymore. Spin-charge separation in 
the usual sensei^ is therefore broken by the SOI. This 
effect can be traced back to a term in the SO Hamil- 
tonian diagonal in sublattice space (see Eq. [1]), which 
was previously overlooked, (ii) We discuss in detail the 
spectral function, a quantity that can directly be probed 
experimentally for SWNTs.— We show how the mixing 
of spin and charge modes due to SOI affects its ana- 
lytic structure and modifies it from the established spin- 



ful LL bchavior-i^i2i. The predicted deviations are small 
but should be observable, (iii) The tunneling density of 
states, and hence most typical quantum transport ob- 
servables, is only weakly affected by the SOI. This may 
explain why the SOI in SWNTs has long been overlooked, 
(iv) We shall clarify the similarities and the differences 
of the present SWNT theory to the LL description of ID 
interacting semiconductor wires with Rashba SOI<^"— 

To start, let us address the bandstructure of a nomi- 
nally metallic (n, m) SWNT, where 2n + to e 3Z. The 



chiral angle^^ is 



= taii-^[V3m/{2n 
0.0391v^ 



to)], and the 

tube radius is i?[nm] ~ O.OSOlVn^ + nra + to^. We em- 
ploy the effective SO Hamiltonian for tt electrons de- 
rived in Ref. [ll| in the k ■ p scheme. This model is in 
semi-quantitative accordance with available experimental 
Coulomb blockade spectroscopy data,— and summarizes 
earlier theoretical work. In particular, it includes the re- 
cently discovered "diagonal" contribution Eso, which is 
of crucial importance in our analysis (see below). Within 
this framework, the single-particle Hamiltonian Ho{k) for 
wavevector k — (fc, fcj_) relative to the respective K point 
is a 2 X 2 matrix in sublattice space corresponding to the 
two basis atoms of the honeycomb lattice. This sepa- 
rately applies to both K points a = ± and both spin di- 
rections (T = ±, where the spin quantization axis is along 
the tube axis. To leading order in the SOI, the spin label 
a is still a good quantum numberj^i^i^ Periodic bound- 
ary conditions around the SWNT circumference imply a 
quantization of transverse momentum, kj^R = ng G Z. 
We assume a Fermi energy Ep > but sufficiently small 
to justify that only the no = band has to be retained. 
All other bands are then separated by an energy gap 
« hvp/R « 1 eV, where wj. « 8 x 10^ m/s. Then, Ho{k) 
is given byii 



aaEso 

-ahvp[(f>± — i{k 



^011 )] 



i)± + i{k - 
aaEso 



(1) 

This form neglects trigonal warping corrections,— which 
cause only tiny changes in the low-energy physics but 
would complicate our analysis substantially. Using the 



2 



parameter estimates of Ref. [ijj, the diagonal term is 

0.135 cos(36') 



i;so[nieV] 



i?[nml 



(2) 



Writing 4>± = 4>±^so + </'_L,cur, the SOI corresponds to a 
spin-dependent shift of the transverse momentum,— 



2.7 X 10- 



(3) 



i?[nm] 

while curvature effectsiiii^ give 

, 0.011 cos(36') 
0±,curnm ~ —— -2—, 

, , 0.045 sin(36l) 

(/[[nmjj^ 

We remark that the Hamiltonian ([T]) contains the two 
leading effects of curvature-induced SOI, namely the di- 
agonal contribution Ego and the Rashba-type SOI en- 
coded by 0±,so- Subleading terms, e.g., the "intrinsic" 
SOI)^ are much smaller and not taken into account here. 
The dispersion relation obtained from Eq. ([ij is 



E^^-''\k)^aaEso±fiVFyJc^i + {k + a(^\\)^ (4) 
where the Kramers degeneracy is reflected in E^''^'^ (k) ~ 



Since Ep > 0, only the conduction bands 

(positive sign) are kept, and the Fermi momenta ki-aa 
for right- and left-movers (r = R/L = ±) follow from 

e'^'''^ (^kiaVj = Ep, kial ~ r{EF - aaEso) /fiVF - a(l)\\. 

We linearize the dispersion relation around the Fermi 
points, always assuming that -E_f is sufficiently far away 
from the band bottom. The ID Fermi velocities Va^a = 

fi~^9fc£'^"''^^ (jt = fc+^Qcr) take only two different values, 

va = 'I'-.t ~ ''^+,1 aiT^d vb = ~ ''^-.i- We men- 

tion in passing that R/L movers have pairwise identi- 
cal velocities only in the absence of trigonal warping 
and orbital magnetic fields^^ or transverse fields )2£ as 
assumed here. It is convenient to introduce the mean 
velocity v = {va + wb)/2 and the dimensionless differ- 
ence 6 = {vA — vb)/{2v). After some algebra, Eq. 
together with the parameter estimates above yields 



V 

VF 



0.01(i?[nm])2 + 17cos2(36l) 
{EF[iney]y (i?[nm])4 
0.83cos(36i) 



(£;F[meV])2 (i?[nm])3 



(5) 



The renormalization of v away from vf goes always 
downwards, but the quantitative shift is small. The 
asymmetry parameter S effectively parametrizes the SOI 
strength and is more important in what follows. For 
fixed Ep and R, it is maximal for 9 = (zig-zag tube) 
and vanishes for 6 = t:/6 (armchair tube). Moreover, 



S increases for smaller tube radius, but the continuum 
description underlying our approach eventually breaks 
down for R < 0.4 nm. Since Ep should at the same time 
be sufficiently far above the band bottom in Eq. in 
practice this leads to rather small values, S < 0.05. This 
is a rather conservative estimate, though, based on the 
parameter values of Ref. [ll| and larger values could be 
obtained if one uses different estimates. Nonetheless, we 
show below that observable consequences do arise. 

The theory is then equivalently formulated using 
Abelian bosonizationji^ which allows for the nonpertur- 
bative inclusion of interactions. We employ the boson 
fields (paix) with a = c+,c—,s+,s—, representing the 
total and relative charge and spin density modeSfii and 
their conjugate momentum fields Ila{x) = —dxda, where 
9a arc the dual fields. Those fields arc conveniently com- 
bined into the vectors $/(a;) = (0C+, 6'c+, 4>s-: d.s-)'^ and 
^o{x) = {<j>c-,Sc-,4's+,0s+)'^ ■ The important electron- 
electron forward scattering^ effects are parametrized 
by the standard LL parameter K = Kc+, where 
K = 1 for noninteracting electrons but K sa 0.2 . . . 0.4 
for SWNTs deposited on insulating substrates (or for 
suspended SWNTs) due to the long-ranged Coulomb 
interactioniii"— The low-energy Hamiltonian of a spin- 
orbit-coupled interacting metallic SWNT then reads 



H 



( h[K) 
^ h{l) 



with the A'-dependent matrix 



h{K) 



(6) 




The above representation shows that SOI [6 ^ 0) breaks 
spin SU(2) symmetry. Notably, the modes $/ and $o 
decouple, and interactions [K ^ 1) only affect the $/ 
sector. In each sector, the Hamiltonian is then formally 
identical to the one for a semiconductor wire with Rashba 
SOI in the absence of backscattering<^ We consider a 
very long SWNT and ignore finite-length effects, i.e. the 
zero modes contributions to the Hamiltonian ([H]). 

Equation ([6]) can be diagonalized by the linear 
transformation^ <!>/ = Vi^a and $0 = Vb<i>6, with the 
4x4 matrix 



Vi = 



where 



/ COS77 

cos ?7 —y sin 77 

y sin ?/ cos rj 

0^0 cos?7 / 

25y 



(7) 



= ^{l + K-^)/2, tan(277) 



y 



(8) 



Vq is as in Eq. ([7]) with if = 1, i.e., y = 1 and ?/ = 7r/4. 
In terms of the new vectors $p = (</>+, p, ^+.p, <i>-.pi ^-,p)^ 



<jp and 9jp for each set 



with mutuahy dual boson fields 
(j = ±,p = a/h), the diagonalized Hamiltonian is seen 
to describe a four-channel Luttinger liquid, 



H = 



dx 



1 



K 



3P 



(9) 

The interacting sector corresponds to p = a, where the 
effective LL parameters K±m and the plasmon velocities 
v±,a are 



K±M 

V±.a 
V 

A 



3 + X-2± A 
37^-2 + 1± A' 



= vV+12±a72, 



(10) 



with y in Eq. ([5]). For p = b, the nonintcracting values 
apply, ii'i^b = 1 and v±^i, = v{l ± S). Note that the 
above expressions recover the LL theory for S = Oj^^ 
where Vjp = vp/Kjp with = 1 except for K+^a = K. 

Within the framework of the LL Hamiltonian us- 
ing the bosonized form of the electron field operator-i^^ii 
^ raa{x,t) and the transformation ([7]), it is possible to 
obtain exact results for all observables of interest. In 
particular, arbitrary correlation functions of exponentials 
of the boson fields can be calculated. As an important 
application, we discuss here the spectral function for an 
r = R/ L moving electron with spin a near the K point 
a = ±, which is defined as 



1 



Im G^'=„V('7,w), 



(11) 



with the Fourier transform of the single-particle retarded 
Green's function (0 is the Heaviside function), 

Glt.M = -iQ{t) [(*,„,(a;,t)vl/t^,(0,0)) +C.C.] , 

and the momentum q is measured with respect to the 
relative Fermi momentum fcrao-- 

After some algebra, Eq. ([TT|) follows in closed form, 
which we specify in the zcro-tcmperature limit now. 
With the short-distance cutoff (lattice spacing) ao ~ 
0.246 nm, we find 



tT(g,w) oc 

nn 

■ hP A'=± 



dx 



.Vnpt + 



dt 



p.rx\ 



where the exponents for = ± and p, 
[see also Eq. ([8])] 



(12) 

-I- {x,t) {-x,-t) 
± are given by 



1 r 

16 



cos(7?) (^K^l^ - pK^ l 
iv) {y^K]'a - p.y' 



(13) 



aa] sm 



-1/2 
J, a 



A(co) 
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FIG. 1: (Color online) Spectral function (|14|l for a right- 
mover in an interacting SWNT with LL parameter K = 0.4 
and SOI parameter 5 = 0.05, shown in arbitrary units as 
function of uj for given wavevector g > 0. The black solid 
curve is for aa = 4-1, while the red dashed curve is for ao — 
— . Note that ARacr{q,i^) ~ for —v-.aq < co < vq. Right 
inset: Magnified view around lj ~ V-,aq- Left inset: Same as 
main panel but without SOI (5 = 0). Shifts of the positions of 
the singularities due to the shifts of Fermi momenta are not 
included in the figure since each spectral function Araa{q,ui) 
is evaluated at momentum q relative to the respective Fermi 
momentum. 



The remaining Fourier integrals are difficult to perform. 
We here follow Ref. [H and focus on the analytic structure 
of the spectral function, which can be obtained by the 
power counting technique and Jordan's lemma. Up to an 
overall prefactor, the spectral function exhibits power- 
law singularities close to the lines u = ±Vjpq. These 
singularities are captured by the approximate form 



Arc^{q,io) « m Icj + ^rwj- agr'""' ^ I (14) 

I ^-1 , f\ ir*°"'-3/2 
X \u! ~ r[l + aad)vq\ 

X [Q{iu — rvq) + 0(— cj — rv^al)] 
where v = min[w_.a, (1 + aa5)v\ and 



(aa) 
jPfi 



(15) 



Wc stress that Eq. is asymptotically exact: it has the 
same analytic structure and the same exponents of the 
power laws at the singular lines uj = ±Vjpq as the exact 
spectral function. Away from the singularities, however, 
it only serves illustrative purposes. 

The spectral function is depicted in the main panel 
of Fig. [T] for fixed wavevector g > as a function of fre- 
quency w, taking K = 0.4 and i5 = 0.05. Compared to 
the well-known spectral function in the absence of SOI 
{S = 0), see left inset of Fig.[T]and Refs. [l6ll2l1 . additional 



4 



structure can be observed for d 0. First, the singular 
feature around uj = V-^aQ splits into two different power- 
law singularities when S ^ 0, sec the right inset of Fig. [T] 
for a magnified view. For large q, the corresponding fre- 
quency differences are in the meV regime and can be 
resolved even for the rather small S expected here. Sec- 
ond, for — 7j_|-.a'? < w < —v^^aQ, the spectral function is 
finite (albeit small) when S 0. Note that for S = 0, the 
respective velocities are f+^a = vp/K and V-^a ~ vp, 
implying a large frequency window where this effect may 
take place. These predictions for the spectral function 
could be detected by photoemission spectroscopy. 

Many standard quantum transport properties, how- 
ever, will hardly show an effect due to the SOI, which 
may explain why effects of SOI in SWNTs have been so 
long overlooked. For instance, the tunneling density of 
states averaged over (r, a, a) exhibits power-law scaling 
with u for low frequencies, v{u]) cx w'*'"^. The exponent 
7 is the smaller of the quantities F*^*) in Eq. ([TS]) . This 
exponent is analytic in 5, and the smallness of 5 then 
implies that the tunneling density of states in SWNTs 
will be very close to the one in the absence of SOI. 
Let us also briefly comment on the relation of our re- 



sults to the LL theory for semiconductor quantum wires 
with Rashba S0I,^"-2^ The "interacting" sector p = a 
in Eq. ^ coincides with the semiconductor theory when 
electron-electron backscattering can be neglected. The 
additional presence of the "noninteracting" sector p = b, 
however, causes additional structure in the spectral func- 
tion. Moreover, while backscattering in semiconductor 
wires is likely an irrelevant perturbation in the renormal- 
ization group sense^^ it nonetheless causes a renormal- 
ization of the LL parameters and the plasmon velocities. 
Such renormalization effects are negligible in SWNTs. 

To conclude, we have studied SOI effects on the effec- 
tive low-energy theory of interacting metallic SWNTs. 
We have shown that a four-channel Luttinger liquid the- 
ory remains applicable, but compared to the previous for- 
mulation without SOI,-i^ all four channels arc now charac- 
terized by different Luttinger liquid parameters and plas- 
mon velocities, reflecting the broken spin SU{2) symme- 
try. The coupling of spin and charge modes leads then to 
observable modifications in the spectral function, which 
provide an experimental signature of SOI. This work was 
supported by the SFB TR 12 of the DFG. 
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